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1  Introduction 

In  the  past  term  we  have  continued  our  effort  towards  understanding  amplification  in 
single  as  well  as  double  wavelength  pumped  parametric  amplifiers  when  these  are 
operated  in  depletion.  Two  applications  that  we  have  in  mind  are  wavelength 
conversion  and  signal  regeneration.  In  addition,  we  have  also  looked  into  the 
possibility  of  using  the  fourth  order  derivative  of  the  propagation  constant  to 
obtain  phase  matching  at  large  wavelength  shifts  from  the  wavelength  of  the  pump. 

Our  effort  has  been  mainly  theoretical  however,  we  are  in  the  process  of  designing 
and  putting  together  a  new  phase  sensitive  amplifier.  The  following  sections  provide  a 
status  report  on  our  results.  In  Section  2  results  on  investigations  of  fluctuations  of  the 
group  velocity  dispersion  on  the  gain  spectrum  of  a  parametric  fiber  amplifier  are 
described.  This  effort  has  resulted  in  a  publication  at  the  forthcoming  OS  A  Topical 
meeting  on  nonlinear  photonics.  The  manuscript  is  included  in  enclosure  1.  In  Section 
3  theoretical  results  on  depletion  and  nonlinear  Stokes  analysis  is  described.  We  are 
currently  putting  these  results  together  in  two  publications  that  we  expect  to  submit  to 
Optics  Express  within  the  coming  months.  Both  papers  at  their  current  state  are 
included  to  this  report  as  enclosures  2  and  3.  Finally,  Section  4,  describes  our  current 
effort. 

Section  5  concludes  this  status  report. 


2  Influence  of  variations  of  group  velocity  dispersion  (enclosure  1) 

The  spectral  shape  of  the  gain  of  a  parametric  process  in  an  optical  fiber  depends  on 
the  group  velocity  dispersion  of  the  optical  fiber,  the  initial  pump  and  signal  powers, 
and  the  wavelength  of  the  pump.  These  parameters  can  be  chosen  to  create  a  balance 
between  the  second  and  fourth  order  group  velocity  dispersion  terms,  such  that  the 
gain  spectrum  consists  of  two  gain  regions  on  each  side  of  the  wavelength  of  the 
pump,  see  Fig  1,  which  show  the  gain  spectrum  resulting  from  a  350  m  long  fiber 
pumped  by  1  Watt.  Three  cases  are  considered  corresponding  to  three  different  pump 
wavelengths.  The  second  gain  region,  i.e.  the  gain  region  furthest  away  from  the 
pump,  is  of  interest  for  broad  wavelength  conversion  or  when  working  outside  the 
conventional  C  and  F  bands.  For  narrow  peaks  (<  1  nm)  it  also  enables  applications 
such  as  a  narrow-band  filter  with  gain  or  as  a  gain  medium  in  fiber  lasers  with 
tuneable  wavelength.  On  the  other  hand  if  the  gain  spectrum  is  broad  (>  10  nm),  it 
could  be  used  for  amplification  in  WDM  systems,  while  still  avoiding  four  wave 
mixing  between  the  channels. 

The  numerical  simulations  shown  in  Fig.l,  illustrate  how  random  fluctuations, 
according  to  a  normal  distribution,  of  the  group  velocity  dispersion  along  the  fiber 
strongly  influence  both  gain  peaks  in  the  gain  spectrum  and  especially  the  gain  peak 
occurring  furthest  away  from  the  pump  wavelength.  The  considered  fluctuations  in 
the  group  velocity  dispersion,  are  similar  to  those  existing  in  state  of  the  art  highly 
nonlinear  optical  fibers.  From  the  figure  it  is  obvious  that  the  spectral  gain  profile  is 
strongly  modified  and  the  second  gain  peak  disappears. 
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Further  details/discussions  on  this  topic  are  found  in  enclosure  1,  whieh  is  a  copy  of 
the  manuscript  accepted  for  presentation  at  the  fortheoming  OS  A  Topical  Meeting  on 
Nonlinear  Photonics,  to  be  held  in  Karlsruhe,  Germany,  June  21-24,  2010.. 


Figure  1:  Gain  versus  signal  wavelength  for  three  pump  wavelengths.  The  launched  signal  and  idler 
power  is  0.1  mW.  The  lines  reaching  25  dB  of  gain  are  with  no  variations  in  group  velocity  dispersion, 
whereas  the  other  lines  are  obtained  with  realistic  variations  in  group  velocity  dispersion. 


3  Other  results 

In  addition  to  the  work  on  impact  do  to  fluctuations  in  group  velocity  dispersion  we 
have  continued  our  effort  on  theoretical  modelling  of  the  parametric  fiber  amplifiers. 
Our  most  recent  results  are  described  below. 

3.1  Semi-analytical  model  (enclosure  2): 

One  important  application  of  the  parametric  amplifier  is  when  it  is  operated  in 
depletion,  since  this  has  use  in  communication  systems  as  regenerating  amplifiers  or 
as  wavelength  converters. 

We  have  derived  a  semi-analytic  model  of  the  dual  pumped  parametric  amplifier  that 
includes  pump  depletion.  The  analytic  solution  is  found  by  using  Jacobian  Elliptic 
functions.  These  functions  are  not  easy  to  use  however,  with  a  few  but  reasonable 
approximations  we  have  obtained  some  much  simpler  and  useful  expressions.  We 
have  used  this  model  to  show  that  wide  and  fiat  gain  is  achievable  with  the  dual 
pumped  FOPA.  In  addition,  the  theoretical  modelling  shows  that  a  phase  mismatch 
which  equals  the  product  of  the  nonlinear  strength  times  the  signal  power  leads  to 
complete  power  conversion. 

Further  details  and  discussions  are  found  in  enclosure  2,  which  is  a  version  of  a 
manuscript  that  we  are  currently  writing,  which  describes  the  semi-analytical  model 
of  the  depleted  parametric  amplifier.  We  plan  on  submitting  the  manuscript  to  Optics 
Express  within  May. 

3.2  Nonlinear  Stokes  analysis  (enclosure  3) 

To  get  further  into  a  detailed  description  of  noise  issues  we  have  studied  the  dynamics 
of  the  four-wave  mixing  (FWM)  equations  using  a  Stokes  like  parametrization  of  the 
Fields.  The  parametrization  is  chosen  so  as  to  use  symmetry  of  the  Gauge  invariance 
of  the  FWM  equations.  This  has  enabled  us  to  reduce  the  dynamics  of  FWM  to  the 
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one-dimensional  interseetion  between  the  phase-plane  of  the  eonserved  Hamiltonian 
and  the  elosed  two-dimensional  iso-power  surfaee  in  a  three-dimensional  phase-spaee. 
In  this  way  both  the  phases  and  amplitudes  of  the  mixing  fields  ean  be  visualized 
simultaneously  giving  a  eomplete  overview  of  the  FWM  dynamics.  Further  details  are 
described  in  enclosure  3,  which  is  a  current  version  of  manuscript  that  we  are 
planning  on  submitting  to  Optics  Express  within  the  coming  months. 

The  analysis  is  advantageous  since  non-degenerate  FWM  as  in  two  pump  optical 
parametric  amplification  is  strongly  phase  dependent.  The  analysis  also  visualize  the 
effect  of  noise  as  noise  will  manifest  itself  as  a  perturbation  of  the  iso-power  surface 
which  might  be  useful  for  analysis  of  the  noise  limits  for  depleted  parametric 
amplifiers. 


4.  Ongoing  activities 

The  activities  described  above  are  all  completed  to  a  level  where  results  are  close  to 
being  submitted  for  publication  or  even  published.  In  addition  to  these  activities  we 
also  have  activities  that  we  are  currently  working  on,  but  where  we  at  present  time  do 
not  have  results  that  may  be  published. 

Most  important  are  our  activities  related  to  systems  demonstrations  of  phase  sensitive 
amplifications.  These  activities  have  been  funded  by  the  Danish  technical  research 
council  but  have  been  seeded  by  this  current  project.  We  have  only  just  started  the 
project  and  have  hired  a  post  doc  and  one  ph.d.  student.  However,  at  this  point  in  time 
we  have  no  results. 

In  addition  to  the  program  on  phase  sensitive  amplifiers  we  have  initiated  a  master 
student  project  entitled:  “Pulse-deformation  in  parametric  amplifier”  Within  the 
project  a  parametric  amplifier  will  be  used  as  the  amplifier  in  a  fiber  ring  laser,  and 
formation  of  short  pulses  will  be  studied.  The  program  has  just  been  initiated  and  at 
this  point  in  time  the  project  has  not  given  any  results. 

Finally,  in  collaboration  with  Colin  Mc.Kinstrie,  Alcatel-Lucent,  NJ,  USA,  we  have 
been  working  on  Self  Seeded  Four  wave  mixing.  This  is  the  process  where  strong 
pumps  interact  during  propagation  and  create  two  additional  self-seeded  pumps.  The 
two  strong  pumps  interact  in  degenerate  parametric  amplification  to  create  two 
sidebands.  When  these  are  created,  all  four  also  interacts  in  a  non-degenerate 
parametric  process.  The  evolution  of  these  four  pumps  is  investigated  with  focus  on 
creating  four  equal  pumps  suited  for  multicasting. 

We  have  performed  a  theoretical  study  of  this  process  and  are  currently  discussing  the 
results. 
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5.  Conclusion 

In  the  term  eovered  by  this  report  we  have  foeussed  partly  on  the  impaet  of 
fluctuations  along  the  length  of  the  fiber  through  which  amplification  is  achieved  in 
the  parametric  amplifier  and  partly  on  theoretical  modeling. 

In  parametric  amplifiers  there  exists  two  gain  peaks  one  relatively  close  to  the  pump, 
tens  of  nanometers  from  the  pump,  and  one  further  away,  up  to  hundreds  of 
nanometers  from  the  pump.  The  second  peak  is  very  interesting  for  the  purpose  of 
wavelength  conversion.  However,  we  have  demonstrated  that  due  to  variations  in  the 
group  velocity  dispersion,  similar  to  those  existing  in  state  of  the  art  highly  nonlinear 
optical  fibers,  the  spectral  gain  profile  is  strongly  modified  and  the  second  gain  peak 
disappears. 

Regarding  the  theoretical  modelling  we  have  shown  a  semi-analytical  model  which  is 
capable  of  predicting  gain  of  single  as  well  as  dual  pumped  parametric  amplifiers 
when  these  are  operated  in  depletion.  From  the  theoretical  modelling  it  is  shown  that  a 
phase  mismatch  which  equals  the  product  of  the  nonlinear  strength  times  the  signal 
power  leads  to  complete  power  conversion. 

5.1  Publications  from  the  project 

L.  S.  Rish0j  and  K.  Rottwitt,  ’’Influence  of  Variations  of  the  GVD  on  Wavelength  Conversion  at 
Second  Gain  Region  of  a  Parametric  Process,”  Accepted  for  presentation  at  Nonlinear  Photonics  2010 

K.  Rottwitt,  J.  R.  Ott,  H.  Steffensen,  S.  Ramachandran,  “Spontaneous  emission  from  saturated 
parametric  amplifiers,”  ICTON’09  Azores,  Portugal,  2009 

C.  Peucheret,  M.  Lorenzen,  J.  Seoane,  D.  Noordegraaf,  C.V.  Nielsen,  L.  Griiner-Nielsen  and  K. 
Rottwitt,  ’’Amplitude  regeneration  of  RZ-DPSK  signals  in  single  pump  fiber  optics  parametric 
amplifiers,”  lEE  photonics  Technol.  Letters,  Vol.  21,  No.  13,  p  872,  2009 

M.  R.  Lorenzen,  D.  Noordegraaf,  C.V.  Nielsen,  O.  Odgaard,  L.  Griiner-Nielsen  and  K.  Rottwitt, 
’’Suppression  of  Brillouin  scattering  in  fibre-optical  parametric  amplifier  by  applying  temperature 
control  and  phase  modulation,”  Electron  Lett.  Vol.  45,  No.  2,  p  125,  2009 

M.  R.  Lorenzen,  D.  Noordegraaf,  C.V.  Nielsen,  O.  Odgaard,  L.  Griiner-Nielsen  and  K.  Rottwitt, 
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C.  Peucheret,  M.  Lorenzen,  J.  Seoane,  D.  Noordegraaf,  C.V.  Nielsen,  L.  Griiner-Nielsen  and  K. 
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6.  Future  work 

People: 

Within  the  period  of  this  report,  Toke  Lund  Hansen  has  started  as  a  post  doc. 
supervised  by  Karsten  Rottwitt.  Toke  will  be  performing  mainly  experimental 
research  on  phase  sensitive  amplification.  In  addition,  a  new  ph.d.  student  Valentina 
Christofioi  has  just  started  studies  toward  a  PhD  degree  on  parametric  fiber 
amplifiers.  In  addition,  activities  within  the  systems  group  at  DTU  Fotonik,  has 
intensified  through  two  new  ph.d  positions  within  parametric  amplifiers,  one  partly 
funded  through  a  research  program  on  phase  sensitive  amplifiers  funded  through  the 
Danish  technical  research  council  and  the  other  funded  by  a  so-called 
internationalization  program. 

Activities: 

As  indicated  in  the  above  DTU  Fotonik  now  has  significant  activities  on  using 
parametric  processes  in  optical  fibers.  This  includes  fundamental  studies  on  noise 
issues,  fibers  for  parametric  devices  and  also  system  aspects  of  parametric  amplifiers 
and  more  specifically  phase  sensitive  amplifiers  and  amplification  of  very  short  pulses 
(few  hundred  femto  second  in  pulse  width).  In  addition,  DTU  Fotonik  also  have 
initiated  collaboration  external  partners  in  relation  to  parametric  amplifiers,  this 
include  collaboration  with  OFS  Fitel  Denmark,  with  Prof.  Stojan  Radio,  University  of 
San  Diego  and  Prof.  Siddharth  Ramachandran  Boston  University. 
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1.  Introduction 

The  shape  of  the  gain  spectrum  depends  on  the  group  velocity  dispersion  (GVD)  of  the  fiber,  the  initial  pump  and 
signal  powers,  and  the  wavelength  of  the  pump.  These  parameters  can  be  chosen  to  create  a  balance  between  the 
second  and  fourth  order  dispersion  terms,  such  that  the  gain  spectrum  consists  of  two  gain  regions  on  each  side  of  the 
wavelength  of  the  pump.  Through  simulations  this  paper  focuses  on  the  gain  peak  occurring  furthest  away  from  the 
pump  wavelength,  referred  to  as  the  second  gain  region.  This  gain  region  is  of  interest  for  broad  wavelength  conversion 
or  when  working  outside  the  conventional  C  and  L  bands.  For  narrow  peaks  (<  1  nm)  it  also  enables  applications  such 
as  a  narrow-band  filter  with  gain  or  as  a  gain  medium  in  fiber  lasers  with  tunable  wavelength.  On  the  other  hand  if  the 
gain  spectrum  is  broad  (>  10  nm),  it  could  be  used  for  amplification  in  WDM  systems  and  hereby  avoiding  FWM 
phase  matching  between  the  channels.  Unlike  previous  publications  within  this  field,  this  model  also  accounts  for 
pump  depletion  and  fiber  loss.  The  former  of  these  have  previously  been  neglected,  by  arguing  that  pump  depletion 
rarely  occurs  in  experiments.  However,  it  is  important  when  trying  to  obtain  large  power  transfer  during  wavelength 
conversion.  Also,  recently  depleted  parametric  amplifiers  have  been  used  for  signal  regeneration  [1]. 

2.  Description  of  the  Model 

A  crucial  criterion  for  the  parametric  process  is  to  achieve  phase  matching  between  the  interacting  waves,  i.e.  between 
the  signal,  idler,  and  pump.  In  the  case  of  a  degenerated  pump,  the  parametric  process  depends  on  the  linear  phase 
mismatch  [2].  By  Taylor  expanding  this  around  the  frequency  of  the  pump,  one  yields 


(1) 


From  this  it  is  seen  that  any  variations  of  the  GVD  along  the  fiber  will  influence  the  gain  spectrum  in  an  unknown 
manner.  The  fiber  parameters  used  for  the  simulations  are  for  a  highly  nonlinear  dispersion  shifted  fiber  (HNLF).  This 


fiber  has  a  nonlinear  coefficient  of  7  =  11.5  W  ^  and  a  loss  of  a  =  0.74  dB/km,  these  are  assumed  to  be 
independent  of  wavelength.  Furthermore,  the  fiber  has  a  zero  dispersion  wavelength  of  Aq  =  1559.5  nm  and  a  slope 


at  this  wavelength  equal  to  S'  =  0.015  ps/ (nm^km).  For  this  HNLF  (3i^p  >  0,  thus  by  ensuring  that  ^2,p  <  0,  it  is 
possible  to  obtain  gain  by  balancing  the  two  first  terms  of  Eq.  (1).  Simulations  have  shown  that  if  no  variation  of  the 
GVD  is  included  it  is  possible  to  obtain  wavelength  conversion  over  more  than  300  nm  by  using  merely  1  W  of  pump 
power,  by  selecting  the  wavelength  of  the  pump  accordingly. 

Variations  of  the  GVD  in  fibers  are  caused  by  changes  in  the  core  radius  and  the  refractive  index  along  the 
fiber.  The  latter  contribution  has  been  neglected  in  the  model,  furthermore  the  effect  these  core  changes  has  on  the 
effective  area  and  thus  the  nonlinear  coefficient  has  also  been  neglected.  The  relative  change  in  the  core  diameter  is 
for  recent  fabricated  HNLF  reported  to  be  around  ±0.03  %  over  500  m  [3].  Using  a  full  vectorial  mode  solver  it  was 
found  that  a  change  of  0.03  %  in  the  core  diameter  of  the  HNLF  described  above,  leads  to  a  change  in  the  GVD  of 
0.03  ps/(nm  •  km)  at  Aq;  this  coresponds  to  a  change  in  the  zero  dispersion  wavelength  of  2  nm.  Furthermore,  the 
simulations  with  the  mode  solver  verified  that  it  is  only  a  minor  asumption  to  model  the  variations  of  the  GVD  as 
vertical  translations  of  the  entire  GVD  curve,  hence  this  is  utilized  in  the  following  simulations.  The  GVD  along  the 
fiber  is  modeled  by  adding  variations  according  to  a  stochastic  variable  following  a  normal  distribution,  this  is  similar 
to  the  method  used  in  [4,  5].  The  distance  between  these  random  generated  points  is  denoted  L„ar  and  the  GVD 
between  these  points  are  aproximated  to  be  linear.  The  symbol  used  for  the  standard  deviation  is  a. 


3.  Simulation  Results 


It  is  found  that  the  pump  is  fully  depleted  when  the  pump  wavelength  is  Ap  =  1564  nm  and  no  variations  of  the  GVD 
are  included,  if  the  signal  wavelength  is  As  =  1554  nm  and  the  fiber  length  is  Lopt  =  347  m,  hence  working  in  the 
first  gain  region.  This  situation  is  referred  to  as  ideal  and  has  been  used  for  normalization.  Since  a  stochastic  variable  is 
included,  the  simulations  have  been  carried  out  1000  times.  It  was  seen  that  the  distribution  for  any  parameters  follows 
an  exponential  behavior,  which  peaks  at  a  normalized  idler  power  of  one.  However,  the  spread  becomes  narrower  as 
Lyar  and  cr  decrease,  these  results  have  not  been  included  here.  In  Fig.  1  it  is  seen  that  the  minimum  power  of  the  idler 
increases  as  Lyar  and  a  decrease.  This  indicates  that  even  in  the  worst  case  scenario  a  significant  conversion  efficiency 
is  obtained,  if  the  variations  are  on  the  order  of  or  shorter  than  a  few  meters. 


Fig.  1.  The  paiameters  for  the  simulation  are  As  =  1554  nm  and  Ap  = 
1564  nm.  The  minimum  normalized  idler  power  after  the  optimal  length 
of  propagation.  The  bottom  10%  of  the  simulations  have  been  removed 
in  order  not  to  obscure  the  tendensies.  The  legend  indicates  the  standard 
deviation  of  the  normal  distribution  in  units  of  ps/(nm  ■  km). 


Fig.  2.  The  gain  spectrum  for  different  wavelengths  and  with  a  fiber 
length  of  L  =  350  m.  The  parameters  are  Pp^i  =  1  W  and  Ps^i  = 
0.1  mW.  The  lines  reaching  gain  over  25  dB  are  with  no  variations  of 
the  GVD,  whereas  the  remaining  lines  are  for  variations  where  cr  = 
0.01  ps/(nm  •  km)  and  Lyar  =  1  m. 


Next,  a  similar  analysis  is  performed  for  the  gain  region  that  occurs  further  away  from  the  wavelength  of  the  pump, 
also  referred  to  as  the  second  gain  region.  In  these  simulations  As  and  L  was  chosen  such  that  full  depletion  of  the  pump 
is  obtained,  if  no  variations  of  the  GVD  are  included.  It  was  seen  that  hardly  any  idler  power  is  generated.  In  general  the 
tendency  is  that  as  Lyar  decreases  or  cr  increases  the  distribution  becomes  narrower  and  closer  to  zero.  Also  the  spread 
of  the  idler  powers  becomes  larger  as  Ap  becomes  smaller,  these  results  are  not  included  here.  Fig.  2  shows  the  gain 
spectra  for  three  different  wavelengths  of  the  pump,  both  with  no  variations  of  the  GVD  and  with  six  different  random 
variations,  the  initial  signal  power  is  set  to  Ps,o  =  0.1  mW.  It  is  seen  that  the  gain  spectrum  is  seriously  degraded  by 
the  introduction  of  variations,  especially  it  is  noticeable  that  the  second  gain  region  is  basically  non-exciting  for  all  of 
the  three  wavelengths  of  the  pump.  However,  a  certain  amount  of  gain  is  still  obtained  for  wavelengths  close  to  the 
wavelength  of  the  pump.  This  is  not  obvious  by  observing  the  Taylor  expansion  of  the  phase  matching  condition,  since 
as  Ap  gets  closer  to  Aq  the  relative  change  in  j32,p  due  to  the  same  sized  variations  gets  larger,  which  would  increase 
the  impact  of  the  variations.  On  the  other  hand  the  gain  peak  is  moved  such  that  \ujp  —  Wsl  becomes  smaller  and  also 
the  gain  peak  becomes  wider,  which  both  reduce  the  impact  from  the  variations.  A  solution  would  be  to  increase  the 
pump  power  and  thereby  reduce  the  needed  fiber  length,  and  thus  in  effect  reduce  the  impact  of  the  variations  of  the 
GVD. 


4.  Conclusion 

It  was  shown  that  introducing  variations  of  the  GVD  of  a  magnitude  similar  to  those  expected  in  modern  HNLF,  the 

parametric  process  was  influenced  to  an  extend  such  that  the  gain  peak,  caused  by  the  balancing  of  the  terms  involving 

I32.P  and  /34  p  in  the  phase  matching  condition,  would  cease  to  exist. 
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Abstract:  This  paper  solves  the  four  coupled  equations  describing 

non-degenerate  FWM,  with  the  focus  on  amplifying  a  signal  in  a  fiber 
optical  parametric  amplifier  (FOPA).  Based  on  the  full  analytic  solution, 
a  simple  approximate  solution  describing  the  gain  is  developed.  The 
advantage  of  this  new  approximation  is  that  it  includes  the  depletion  of  the 
pumps,  which  is  lacking  in  the  usual  quasi-linearized  approximation.  Using 
this  approximate  solution  we  demonstrate  the  possibility  of  achieving  a  flat 
gain  over  87  nm  centered  at  1560.5  nm. 
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1.  Introduction 

Parametric  amplifiers  have  shown  promising  features  in  a  wide  variety  of  applications  beside 
just  amplification  of  signals.  It  has  potential  as  optical  regenerators  [1,  2,  3]  based  on  either  one 
or  two  pumps,  to  suppress  level  fluctuation  [4],  and  also  as  noise  limiters,  for  instance  to  limit 
amplified  spontaneous  emission  (ASE)  noise  or  to  reduce  nonlinear  phase  noise  [5]. 

Dual  pumped  Parametric  amplifiers  show  potential  application  as  amplifiers  with  a  very  wide 
bandwidth  of  operation  and  a  flat  gain  throughout  the  bandwidth  window  of  operation  [6]. 
However,  it  is  more  complicated  to  predict  their  performance  analytically  since  four  coupled 
wave  equations  needs  to  be  solved.  In  this  work  we  demonstrate  the  existing  analytical  solutions 
using  elliptic  Jacobian  functions  [7,  8],  and  based  on  these  we  show  that  simple  approximations 
exist  that  enable  predictions  of  depletion  performance  of  double  wavelength  pumped  parametric 
amplifiers. 

Previously  Kylemark  et  al.  [9]  have  derived  a  useful  approximation  for  a  single  pump  FOPA 
operated  at  the  total  conversion  phase  mismatch.  This  paper  uses  a  similar  approach  on  the  dual 
pumped  FOPA  and  extends  it  to  the  full  range  of  phase  mismatch.  We  use  the  approximation 
on  a  specific  example  where  it  is  shown  that  a  flat  and  wide  gain  is  achieved  with  pumps 
placed  symmetric  around  the  zero  dispersion  wavelength  (ZDW).  A  flat  gain  over  25  nm  have 
previously  been  experimentally  realized,  [10],  but  in  this  paper  we  predict  a  flat  broadband  gain 
over  87  nm. 

2.  Theory 

In  general  we  consider  an  electric  field  consisting  of  four  CW  waves,  at  frequencies  (Oi  through 
CO4.  In  non-denegerate  FWM,  the  four  distinct  waves  interact  with  each  other  under  the  condi¬ 
tion  that  Oi  +  0)4  =  ©2  -b  ©3 .  Two  of  these,  number  2  and  3,  are  in  this  paper  used  as  pumps, 
while  1  is  a  signal,  also  denoted  s,  that  is  to  be  amplified,  and  4  is  the  idler,  i,  that  arise  when 
the  signal  is  amplified.  The  envelope  of  the  field  is,  for  parallel  pumps  and  signal,  governed  by 
[11,  12] 


dA 

=  ipi^i  +  )7{  [1^1 1^  +  2  (IA2P  +  IA3 p  +  IA4P)]  Al  +  2A2A3A4}  (la) 

=  /^2^2  +  iy {  [1^2!^  +2  (|Ai  p  +  IA3P  +  IA4P)]  A2  -I-2A1A3A4}  (lb) 

+  /7I  [|A3|^  -1-2  (|Ai  p  +  IA2P  -b  IA4P)]  A3  -b2AiA2A4}  (Ic) 

dAi 

=  /^4A4-b/7{  [|A4|^ -b2  (|Ai  p -b  |A2p-b  IA3  p)]  A4 -b  2A]‘A2A3  }  ,  (Id) 

where  7  is  the  nonlinear  strength  and  j3a  is  the  wave  number,  index  a  may  be  either  1,  2,  3  or  4. 
In  the  case  of  perpendicular  pumps  and  sidebands,  the  factor  of  2  in  the  coupling  term  should 
be  changed  to  a  1  [13,  14].  The  most  common  method  of  solving  Eq.  (1)  is  by  assuming  that 

the  pump  powers  remain  much  greater  than  the  signal  and  idler  power  at  all  times  and  that  the 


power  loss  of  the  pumps  is  negligible  [15].  This  will  quasi-linearize  the  differential  equations 
and  for  pumps  with  identical  power,  Pp,  result  in 
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where  =  {lyPpf-  —  (k'/2)^  and  K  =  27Pp  +  Aj3.  Aj3  is  the  phase  mismatch  given  as  Aj3  = 
j3i  —  ft  —  ft  +  [34,  and  since  Pp  is  the  pump  power  in  each  pump  and  these  are  assumed  identical, 
the  expression  resembles  the  expression  known  from  single  pumped  parametric  amplifiers  [15]. 
Based  on  this  solution,  the  highest  growing  gain,  G(z)  =  Ps{z)/Ps{0),  is  obtained  with  K  —  0, 
thus  requiring  that  Aj3  =  —2YPp.  Generally,  gain  exist  when  —byPp  <  AP  <  2yPp.  The  problem 
with  this  solution  is  that  it  does  not  account  for  either  pump  depletion  or  nonlinear  detuning, 
which  is  caused  by  the  change  in  relative  phase  as  the  power  is  transferred. 

Instead  of  rewriting  the  complex  amplitude  A  in  Eq.  (1)  into  the  real  variables,  power  P  and 
phase  (j),  as  Aq  =  ^/Pi^exp  (i^a),  the  differential  equations  are  instead  written  as 


-^=Sa4y^J[]Pl3sm{e)  (4) 

^=ft  +  y(2^P^-P„)+2y^p^^^cos(e)  (5) 
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where  0  =  +  04  —  ft  —  03  and  Sa  =  It  is  seen  that  only  the  relative  phase,  0,  and 

not  the  individual  phases,  is  important,  thus  the  phase  equations  can  be  reduced  to  one  single 
differential  equation.  As  seen  from  Eq.  (4),  this  relative  phase  is  important  as  it  shows  to  which 
direction  the  power  flows,  if  the  relative  phase  is  positive,  the  power  is  transferred  from  the 
pumps  to  the  signal  and  idler,  if  the  relative  phase  is  negative,  the  opposite  happens.  As  power 
is  transferred  from  the  pumps  to  the  signal  and  idler,  nonlinear  detuning  occurs  which  reduces 
the  coupling.  At  some  point  the  relative  phase  can  be  changed  so  much  that  the  power  flow  is 
reversed  and  the  power  begins  to  couple  back  to  the  pumps. 

Because  the  parametric  process  is  symmetric,  the  signal  and  the  idler  have  the  same  increase 
in  power,  while  the  pumps  have  the  same  decrease,  also  written  as  ft  —Pi)—Q  and  ft(P2  — 
P-i)  =  0.  Einally,  the  process  is  power  conserving  thus  Y.aPa{z)  =  Ptotai-  Based  on  this,  all 
powers  are  related  to  each  other,  and  can  thus  be  described  by  only  one  variable  according  to 
the  Manley-Rowe  relations  [16,  17].  This  variable,  F,  is  chosen  to  be  the  power  of  the  idler 
which  is  assumed  zero  at  z  =  0,  i.e. 


Pi=P,+F 
P2=Pi=Pp-F 
P4  =  F, 


(6a) 

(6b) 

(6c) 


where  Pj  is  the  input  power  of  the  signal,  and  Pp  is  the  input  power  for  each  of  the  two  pumps. 
With  this  notation,  the  gain  of  the  signal  is  found  as 
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From  Eq.  (6),  two  governing  differential  equations  can  be  set  up  for  the  power,  F,  and  the  phase, 

0, 


dF 

dz 

dd 

dz 


=Ay{Pp  -  F)  v'F(P,+F)sin(0) 


=27 


iPp-F) 


-2^F{P,+F) 


cos (0) 


+  75  -  47F 


(8) 


(9) 


where 

5  =  ^+2Pp-P„  (10) 
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This  effective  phase  mismatch  also  includes  the  signal  power,  which  in  the  quasi-linearized  case 
is  omitted  due  to  the  assumptions  made.  The  two  differential  equations  have  the  Hamiltonian 

H  =  AY{Pp-F)^F(P,+F)co&{e)  +  Y(5-2F)F.  (11) 


This  has  been  found  by  using  a  qualified  guess  and  knowing  it  has  to  fulfill  Hamilton’s  equations 
as  given  by 

dF  _  dH  dO  _  dH 

lz~~l)e^  lz~l)F' 


The  first  term  on  the  RHS  of  Eq.  (11)  is  due  to  EWM,  while  the  second  contains  the  linear 
and  nonlinear  phase  terms.  Because  the  Hamiltonian  is  independent  of  z,  it  is  constant.  It  is 
therefore  used  to  reduce  the  problem  to  one  variable.  Since  F  is  initially  zero,  the  value  of  the 
Hamiltonian  is  also  zero.  With  this  the  phase,  0,  can  be  determined  as 


cos (0) 


-{5-2F)F 

4{Pp-F)y¥{P^' 


(13) 


Erom  Eq.  (8)  a  potential  equation  may  be  obtained  for  F,  by  squaring  the  equation  and  inserting 
the  expression  for  0,  resulting  in 


=f[l6iPp-F)^FiP,+F)-i5-2F)^F^ 


For  any  given  value  of  the  wave-number  mismatch,  this  can  be  written  as 


=  l2fF{fi-F){f2-F){F-h) 


(14) 


(15) 


where  fj,  (j  G  (1,2,3))  are  the  roots  of  the  polynomial  and  f\>  f2>  F{Q)  >  f^,.  The  ordering 
of  these  roots  are  important  because  it  affects  the  structure  of  the  solution.  The  solution  for  F 
will  be  bounded  by  the  roots  of  the  potential  equation.  The  lower  boundary  is  zero,  and  the 
upper  boundary  is  the  root  /2.  The  value  of  this  root  will  thus  be  the  maximum  value  of  the 
added  power  to  both  the  signal  and  the  idler.  The  differential  equation  can  be  solved  using 
Jacobian  Elliptic  Eunctions  [18],  the  solution  is: 


F{z) 


/s/zsn^ 

'Yz^/3f]_{f2-  h),m 

-/2  +  /3+/2Sn2 

Yz^/'ifi{f2-f?.),m 

(16) 


where  sn  is  a  Jacobian  Elliptic  function  and  the  squared  modulus 

/2(/l-/3) 


= 


(17) 


hifi-hV 

It  is  possible  to  achieve  complete  power  conversion  from  the  pumps  to  the  signal  and  idler. 
In  this  case  the  power  of  the  idler  equals  the  initial  power  of  the  pumps,  thus  F  =  Pp.  When 
inserted  into  the  Hamiltonian,  knowing  that  the  value  of  the  Hamiltonian  is  zero,  one  see  that 
the  requirement  of  complete  power  conversion  is  that  d  —  2Pp,  giving  a  wave-number  mismatch 
of 


Aj3  =  yPs. 


(18) 


This  differs  significantly  from  the  usually  desired  value  of  —2yPp  which  results  in  the  initially 
fastest  growing  signal.  From  Eq.  (2),  the  initial  growth  rate  for  Aj3  =  —2'YPp  is  G(z)  =  1  + 
sinh^  pyPpz],  whereas  the  growth  rate  for  Aj3  =  yPs  is  G(z}  ~  1  +  |sinh^  [v/jyPpz] .  With  Aj3  = 
yPs  the  roots  of  Eq.  (9)  are  f\  =  fz  =  Pp  and  =  — 4Es/3.  In  this  case,  m  =  I,  reducing 
the  Jacobian  Elliptic  Function,  sn,  to  hyperbolic  tangens  functions,  tanh.  The  solution  is  thus 
rewritten  as 


F{z)  = 


4PpP,sinh2  [yz^Pp{3PpFAP,) 


3Pp  +  4P,  +  4P,sinh2  [yz^Ppi3Pp+4P,)  ' 


(19) 


In  the  interval  —6yPp  <  AP  <  2yPp,  in  which  g  is  a  real  parameter,  approximate  solutions  to 
the  roots  of  Eq.  (14)  is  found  by  treating  P^  as  a  small  perturbation  to  the  potential  equation. 
Doing  this  result  in  three  non-zero  roots. 


=  —l6Ps- 


iM 


rb=Pp 


=  Pp- 


lAp 

6  7 

1  Aj3 

2y 


-12R 


p _ 4P. 

A  y 

ip  _  40 

^^P  y 


-4P 


3Pp  +  2^ 


41 

r 


6Ep-f  y 


(20) 


(21) 


(22) 


The  first  root,  is  the  lowest  root,  thus  it  is  equivalent  to  f^,  and  it  is  seen  that  this  root 
is  proportional  to  —Ps-  The  other  two  roots  are  equivalent  to  f\  and  /2,  however  as  r*  is  not 
always  greater  than  rc,  it  will  depend  on  the  values  of  Aj3,  which  root  is  equivalent  to  which  /. 
However,  the  important  thing  is  to  see  that  they  are  both  proportional  to  Pp.  Therefore,  under 
the  assumption  that  Pp  S>  Es,  m  «  1,  the  elliptic  function  sn  can  be  approximated  with  tanh  for 
small  values  of  z.  Thus,  the  solution  can  be  expressed  as 


F{z) 


f  sinh^ 

h-h 

7^3/1  (/2-/3) 

7^3/1  (/2-/3) 

(23) 


For  very  small  values  of  z,  the  second  term  in  the  denominator  can  be  neglected,  reducing  it  to 
an  expression  similar  to  Eq.  (3)  Consequently,  we  may  apply  the  relations: 


y/l(/2-/3)«^ 


/2/3 

/2-/3 


^yPp 


(24) 

(25) 


which  for  the  interval  —^jPp  <  AP  <  yPp,  have  a  relative  error  of  less  than  1  %  for  P*  =  0  dBm, 
and  this  error  decrease  with  decreasing  input  power.  With  this  F  may  be  expressed  as 


P(z) 


Fql{z) 


1  + 


Fql(z) 


=  F 


^sinh^  [gz] 


1  2 


1  + 


2rPp 


sinh2  [gz] 


(26) 


where  =  (lyPp)^  —  K  =  lyPp+Ap,  and 


1 /M 


Psat  =  fi  • 


{‘^-Ps)+Pp.  -6Pp<f<P, 


R<^<  2P„ 


(27) 


This  final  expression  for  P  is  a  simple  expression,  using  the  quasi-linearized  result  to  create 
an  expression  that  includes  the  saturation  of  the  pumps.  Psat  is  the  upper  limit  of  the  power 
that  can  be  coupled  into  the  idler.  From  this  it  is  also  clear  that  using  a  linear  phase  mismatch 
that  equals  the  nonlinear  phase  shift  per  unit  length  of  the  signal,  yPs,  enables  complete  power 
transfer  from  the  pumps  into  the  signal  and  idler,  while  when  using  a  linear  phase  mismatch 
identical  to  the  nonlinear  phase  shift  of  the  pump  i.e.  Aj3  =  —2yPp  only  can  convert  about  67%. 
However,  in  the  latter  case  when  Aj3  =  —2yPp  the  initial  coupling  efficiency  is  higher. 

If  the  two  pumps  have  different  power  called  Pa  and  where  Pa  >  Pt,  then  all  that  is 
necessary  to  change  is  to  replace  [Pp  —F)^  with  {Pa  —  F){Pi,  —  F)  in  Eq.  (14),  and  2Pp  with 
Pa  +Pb  in  Eq-  (10).  When  this  is  done,  the  solution  is  the  same  as  Eq.  (16).  Because  of  the 
different  pump  powers,  the  phase  mismatch  resulting  in  full  power  transfer  from  the  smallest 
pump  is  Aj3  =  yPs  +  y{Pa  —Pt),  whereas  the  fastest  growing  gain  is  when  Aj3  =  —y{Pa  +Pb)- 


3.  Discussion 

How  well  Eq.  (26)  describes  the  evolution  of  the  idler,  compared  to  Eq.  (2)  is  shown  in  Eig.  1 . 
Erom  the  figure  it  is  seen  that  the  new  expression  fits  well  until  the  power  begins  to  couple  back 
into  the  pump.  In  all  cases  the  new  approximate  expression  has  a  better  correlation  with  the 
analytic  solution  than  the  conventional  quasi-linear  result  in  Eq.  (2). 

Eig.  2  shows  the  gain  for  different  input  powers.  The  gain  has  been  evaluated  by  using 
Eqs.  (2),  (23)  and  (26)  and  it  also  shows  that  the  new  approximation,  Eq.  (26),  is  a  better 
approximation  to  Eq.  (23)  than  the  usual  approximation,  Eq.  (2).  The  new  approximation  takes 
the  saturation  into  account,  but  since  it  does  not  include  back  coupling,  it  is  only  useful  until 
back  coupling  becomes  significant.  Before  the  system  is  in  depletion,  the  gain  is  symmetric 
around  JC  =  0  as  predicted  by  Eq.  (2),  but  when  saturation  set  in,  the  gain  profile  becomes 
asymmetric  which  the  new  approximation  also  takes  into  account.  If  a  setup  is  constructed  with 
two  pumps  placed  symmetric  around  the  ZDW,  then  it  is  possible  to  obtain  a  flat  gain  spectrum 
over  a  range  of  wavelengths,  with  very  high  signal  gain.  The  phase  mismatch  is  estimated  by 
Taylor  expanding  all  the  propagation  constants  around  the  ZDW.  Because  of  the  chosen  pump 
symmetry,  all  where  n  is  odd,  will  cancel  out,  and  since  j32  is  zero  at 

the  ZDW,  the  first  term  that  has  influence  is  j34  which  for  this  setup  is  8.8197  x  10-^®  s^/m. 
The  phase  mismatch  is  thus  approximated  as 

^  [(®s  —  ®zd)"^  —  (®p  —  ®zd)"^]  (28) 

When  using  symmetrically  placed  pumps  and  the  value  of  j34  is  positive  [19],  it  is  very  difficult 
to  obtain  Aj3  =  —2yPp,  unless  the  pumps  are  placed  very  far  from  the  ZDW,  however  Aj3  is 
close  to  zero  for  a  large  range,  resulting  in  a  high  possible  power  transfer  in  this  region.  Eig.  3 


(a)  Ai3/7=  -2P„.  (b)  Ali/Y=0.  (c)  Aj3/7=  ft. 


Fig.  1.  Evolution  of  the  power  of  the  idler  with  Pp  =  30  dBm  and  Ps  =  0  dBm.  The  solid  line 
shows  the  full  analytic  solution,  the  dash-dotted  is  the  usual  quasi-linearized  solution  while 
the  dashed  is  the  approximate  solution  found  in  this  paper.  It  is  evident  that  the  setup  in 
(c)  results  in  complete  power  transfer  while  the  usual  desired  setup  (a)  results  in  the  fastest 
growing  idler,  (b)  shows  that  a  zero  wave-number  mismatch  results  in  almost  complete 
power  transfer. 


(a)  P,  =  -20  dBm 


(b)  ft  =  -10  dBm 


(c)  ft  =  0  dBm 


Fig.  2.  Gain  calculated  at  yL  =  3  with  Pp  =  30  dBm.  The  solid  line  is  the  analytically 
calculated  gain,  the  dashed  line  is  the  gain  calculated  with  Eq.  (2)  while  the  dash-dotted 
line  is  calculated  with  the  approximation  found  in  this  paper.  In  (a)  the  pumps  are  not  in 
saturation,  thus  there  is  no  significant  difference  between  them.  However  in  (b)  the  pumps 
begin  to  be  in  saturation,  which  Eq.  (2)  does  not  take  into  account.  In  (c)  the  input  power 
is  so  high  that  the  pumps  have  depleted  and  for  some  values  of  Aj8  / y  the  power  has  started 
to  couple  back  into  the  pumps  resulting  in  a  lower  gain,  only  the  full  analytic  solution 
accounts  for  this  phenomenon. 


shows  gain  profiles  for  a  fiber  with  L  =  250  m,  Azd  =  1560.5  nm,  and  y=  11.5  W  *km  * 
calculated  with  Eq.  (26). 

Because  the  pumps  are  placed  symmetric  around  the  ZDW,  it  is,  as  seen  in  Fig.  3,  possible  to 
achieve  a  high  and  flat  gain  over  a  broad  spectrum.  Even  if  the  symmetry  is  slightly  broken,  the 
gain  remains  reasonably  flat,  and  even  increase  the  spectrum.  If  Eq.  (1)  is  solved  numerically 
and  including  loss  of  0.74  dB/km,  the  same  spectrum  is  achieved,  only  the  gain  is  1.5  to  2  dB 
lower. 


(a)  Symmetrically  placed  pumps 
with  Xp  —  1564  nm  and  1557  nm. 
The  signal  power  is  -20  dBm,  thus 
the  pumps  are  undepleted. 


Wavelength  [nm] 

(b)  Symmetrically  placed  pumps 
with  Xp  —  1568  nm  and  1553  nm.  the 
signal  power  is  0  dBm,  thus  resulting 
in  depleted  pumps. 


Wavelength  [nm] 

(c)  Asymmetrically  placed  pumps 
with  Xp  —  1564  nm  and  1558  nm  and 
a  signal  power  of  -20  dBm. 


Fig.  3.  Gain  profiles  for  a  case  with  30  dBm  pumps.  With  symmetrically  placed  pumps, 
the  gain  is  very  flat  over  a  large  bandwidth  whereas  with  slightly  asymmetrically  placed 
pumps,  the  gain  is  a  little  less  flat,  but  the  spectrum  becomes  wider. 


4.  Conclusion 

We  have  shown  the  exact  solution  of  a  dual  pumped  FOPA  as  well  as  a  simple  yet  very  useful 
approximation  to  the  solution  which  includes  pump  depletion.  With  the  dual  pumped  scheme  it 
is  possible  to  obtain  a  very  flat  gain  over  a  wide  bandwidth.  With  symmetrically  placed  pumps, 
Fig.  3(a),  the  width  of  the  gain,  within  3  dB  of  max,  is  87  nm,  in  the  case  of  asymmetrically 
placed  pumps.  Fig.  3(c),  the  width  is  increased  to  113  nm.  The  bandwidth  may  be  further 
broadened  by  using  slightly  asymmetrically  placed  pumps.  However  this  is  only  obtainable  at 
the  expense  of  a  less  flat  gain. 

We  have  also  shown  that  although  the  gain  is  increasing  more  rapid  as  a  function  of  fiber  length 
when  Aj3  =  —lyPp,  then  this  phase  mismatch  will  not  result  in  a  complete  power  conversion 
as  is  the  case  when  Aj3  =  yPs- 
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We  consider  the  phenomenon  of  four-wave  mixing.  The  evolution  is  analyzed  using  a  variable 
transform,  making  it  possible  to  visualize  the  power  and  phase  information  simultaneously.  A 
method  of  evaluating  the  dynamics  of  the  system  graphically  without  use  of  complicated  calculations 
is  then  proposed.  The  current  version  is  a  draft  and  therefore  contain  most  calculations. 


I.  INTRODUCTION 

The  optical  parametric  process  of  four-wave  mixing  (FWM)  is  currently  used  for  several  optical  communication 
schemes  such  as  amplification,  signal  copying,  phase  conjugation,  sampling  and  regeneration  [1,  2]. 

The  governing  equations  and  their  evolution  given  by  solutions  to  elliptic  integrals  are  at  this  point  well  known  [3,  4], 
but  the  usual  methods  to  evaluate  power  and  phase  can  seem  complicated.  A  recent  analysis  used  the  assumption  of 
strong  pumps  to  visualize  a  specific  scheme  by  the  use  of  a  coordinate  transform  to  Stokes-like  variables  [5]  following 
a  routine  which  has  previously  been  used  to  describe  the  evolution  of  polarization  [6] . 

In  the  current  work  we  propose  an  easy  way  of  visualizing  the  held  amplitudes  and  relative  phase  of  the  full 
nonlinear  problem  simultaneously  as  dynamics  on  a  closed  2  dimensional  surface  in  a  3  dimensional  space.  Using  the 
method  proposed,  the  dynamics  is  found  without  having  to  evaluate  any  pesky  integrals  nor  differential  equations, 
but  simply  plotting  two  implicit  equations  and  and  observe  their  intersection.  A  similar  approach  has  been  used 
previously  for  three  wave  mixing  in  quadratic  nonlinear  media  [7].  This  showed  a  geometric  interpretation  of  the 
phenomenon  and  gave  a  simple  method  of  visualizing  quasi-phase  matching  used  in  parametric  amplihcation  with 
quadratic  materials. 

The  following  analysis  will  start  by  introducing  the  governing  equations  for  FWM  in  section  II.  Furthermore  this 
section  review  the  structure  of  the  Hamiltonian  and  the  use  of  Poisson  bracket  formalism  for  equating  governing 
equations.  Section  III  introduce  the  geometric  representation  of  the  FWM  equations  using  Stokes  vectors  and  the 
Hamiltonian  in  this  new  space  is  found  and  the  governing  equations  for  the  FWM  in  Stokes  space  are  solved.  In 
section  IV  the  geometric  representation  is  used  to  evaluate  the  dynamics  of  some  optical  communication  schemes 
showing  the  advantages  of  the  proposed  method. 


II.  THE  FOUR- WAVE  EQUATIONS  AND  THEIR  SYMMETRIES 

In  the  following  section  the  four-wave  phenomenon  is  shortly  reviewed  along  with  the  Hamiltonian  and  the  use  of 
the  Poisson  bracket  formalism  with  the  Hamiltonian.  The  well  known  governing  equations  for  the  non-degenerate 
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FWM  process  are  [8] 


^  =  ^7K[|Al|2  +  2(|^2|"  +  I^3|'  +  |^4|')]Ai  +  2^7KA2A3^:e-*^^^  (la) 

^  =  ^7K[|A2|"  +  2(|71i|2  +  |2l3|2  +  |A4|2)]7l2  +  2^7K7la*7l4e*^^^  (lb) 

^  +  2(|^i|"  +  |7l2|2  +  |A4|2)]7l3  +  2^7K7llA;7l4e*^^^  (Ic) 

A  A 

=  iM\Ai?  +  2(|71i|2  +  |7l2|2  +  |A3|2)]7l4  +  2^7K7l^7l27l3e-*^^^  (Id) 


These  are  obtained  by  a  slowly  varying  amplitude  approximation  of  Maxwell’s  equations  for  resonant  interaction  of 
four  lightwaves  propagating  in  a  centrosymmetric  material.  Aj  is  the  slowly  varying  complex  held  envelopes  for  the 
j’th  held,  7k  is  the  nonlinear  Kerr  coefRcient  assumed  to  be  frequency  independent,  and  A/3  =  fii  —  132  —  P3  +  Pa  is 
the  wavenumber  mismatch  with  the  wavenumber  of  the  j’th  held  being  Pj  =  P{ujj).  The  numbering  is  either  from 
high  to  low  frequency  or  from  low  to  high.  The  waves  1  and  4  are  thus  the  side  bands  while  2  and  3  denote  the  inner 
bands. 

By  inspection,  the  set  of  Eqns.  (1)  has  the  Hamiltonian 


4  4 

^  ^  ^  +  27k  {AlA2A:iAle-^^^^ 

n—1  m—1 


AiA*AlAie^^^^) , 


(2) 


with  dn,n  =  5  and  dn,m  =  1  when  n  m. 


Making  use  of  the  conservation  of  power  Pt  =  X^n-i  \  An\‘^  and  introducing  the  change  of  variables  ^ 
the  rescaled  rotating  frame 

=  7k  Pt  and 

Aj  =  qj^/Pte-^^^^+^AKPtz 

(3) 

the  four-wave  equations  become 

^  =  iFigi  -  /|9ip9i  -f  2/^29394: 

(4a) 

^  =  *r292  -  i|92p92  +  2/91^394, 

(4b) 

^  =  «r3g3  -  /|93p  93  +  2/9192  94, 

(4c) 

^  =  ir4g4  -  /|94p94  +  2/g*g293, 

d^ 

(4d) 

with  Fj  =  being  the  rescaled  wavenumber. 
The  Hamiltonian  for  these  normalized  helds  is 


H  =  ^  r„|g„p  “  ^  ^  l9nl"‘  +  2  {q*iq2q3ql  +  91929394)  •  (5) 


Thereby  the  equations  have  the  canonical  Hamiltonian  structure 


c% 

d0 


f)  M 


with  the  Poisson  bracket  dehned  as  [9] 


{P,  Q}qc,q^ 


dP  dQ 
dqa  dql 


dP  9g\ 

dql  dqa  J  ’ 


(6) 

(7) 


with  {P,  P}  =  0  and  {P,  Q}  =  — {Q,P},  which  is  the  classical  counterpart  to  the  quantum  mechanical  commutator. 
The  index  9*  denote  the  choice  of  canonical  variables  for  the  Poisson  bracket  and  will  be  omitted  in  the  following. 
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Due  to  the  normalization,  Eq.  (3),  the  Hamiltonian  no  longer  explicitly  depend  on  the  propagation  parameter,  it  is 
a  constant  of  motion.  Furthermore  the  governing  equations  have  the  conserved  quantities 


=  kiP  -  k4l^ 

(8a) 

i^2  =  |g2p-|g3l^ 

(8b) 

i^3  =  |<?lp  +  k2l^ 

(8c) 

i^4  =  |g3p  +  k4p, 

(8d) 

P=l<7l|"  +  k2|"  +  |g3|"  +  |<Z4^ 

(8e) 

also  known  as  the  Manley- Rowe  relations  (MRR)s  [10,  11].  The  physical  significance  of  these  can  most  easily  be 
understood  by  a  quantum  mechanical  interpretation  as  follows.  The  first  two  of  the  relations  signify  that  in  the  inner 
and  outer  bands,  photons  are  created  or  annihilated  in  pairs.  The  third  and  fourth  relation  signify  that  creating 
the  photon  pairs  in  the  inner  or  outer  bands  annihilates  photon  pairs  in  the  outer  or  inner  bands.  Using  the  first 
two  relations  either  of  the  helds  in  the  third  and  fourth  relation  can  be  interchanged  with  its  held  pair  partner,  i.e. 
K2,  —  Ki  =  |(jf4p  -I-  |(72|^  which  is  the  same  as  interchanging  1  and  4.  The  last  of  the  relations  signify  conservation  of 
photons  and  due  to  the  rescaling  of  the  helds  P  =  1.  The  MRRs  are  seen  to  be  a  overfull  set  so  that  only  three  of 
them  are  independent. 

Furthermore  the  Hamiltonian,  and  thus  the  governing  equations  for  the  helds,  is  seen  only  to  depend  on  the  phase 
diherence  Acj)  =  4’2  +  4>3~4’i~4’i^  where  the  phase  is  dehned  by  (^)  =  | gj (^ | .  Therefore  the  eight  coupled  real 
nonlinear  differential  equations  governing  FWM  can  be  reduced  to  four  equations  governing  the  real  amplitudes  and 
one  governing  equation  for  the  phase  difference.  This  can  be  reduced  to  a  single  potential  equation  for  either  of  the 
amplitudes  which  is  solvable  using  elliptic  integrals.  Subsequently  all  the  other  amplitudes  and  the  phase  difference 
can  be  calculated  [3,  4].  In  the  following  the  problem  will  be  solved  in  another  way  giving  a  geometric  interpretation 
of  all  the  powers  and  the  phase  difference  simultaneously.  The  current  solution  method  also  include  the  common  one 
as  a  special  case. 

As  stated  by  Noether’s  Theorem  the  four  invariant  MRRs  correspond  to  four  symmetries  of  the  governing  equations, 
Eqs.  (4).  These  symmetries  are  the  Gauge  invariances 


(91, 92,  g3,  <74)  - 

(9ie"*'^G92,93,94e*'^"), 

(9a) 

(91, 92,  g3, 94)  - 

(91, 926"*'^",  936**^",  94), 

(9b) 

(91,92,93,94)  - 

(916"*'^",  926“*'^",  93,  94), 

(9c) 

(91,92,93,94)  - 

(9l,92,93e“*‘^^94e"*‘^‘‘). 

(9d) 

This  is  seen  by  letting  either  of  the  MRRs  act  as  the  Hamiltonian  in  Eq.  (6)  and  determining  the  effect  of  the 

evolution,  i.e.  calculating  ^  =  {qj,Kk}.  E.g.  for  Ki  this  would  result  in 

‘^91  T  1  dATi  . 

—  =  i{qi,Ki}  =  =  igi,  (10a) 

=  (10b) 

^  =  ife,A'i}  =  >^=0,  (10c) 

^  =i{q4,Ki}  =  =iq4,  (lOd) 

aq^ 

i.e.  an  equal,  but  opposite,  phase  shift  of  gi  and  g4  which  corresponds  to  the  Gauge  invariance  Eq.  (9a). 

It  is  easy  to  show  that  the  four- wave  equations,  Eqs.  (4),  are  invariant  under  the  Gauge  transforms,  Eqs.  (9). 


III.  REDUCTION  TO  FOUR- WAVE  SURFACES 

In  this  section  the  dynamics  of  the  four-wave  equations  will  be  reduced  to  dynamics  on  a  closed  surface  in  a  3 
dimensional  space.  A  simple  example  of  previous  use  of  such  reduction  is  that  of  projection  of  polarization  onto 
the  Poincare  sphere.  In  that  case  the  reduction  of  dimension  is  performed  by  using  a  set  of  invariant  coordinates 
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with  respect  to  the  phase  symmetry  of  the  polarization,  the  Stokes  parameters,  that  project  the  dynamics  onto  the 
Pointcare  sphere.  Similar  in  our  case  a  set  of  Stokes-like  parameters  which  are  invariant  under  the  same  Gauge 
transforms  must  be  chosen.  Such  a  set  of  real  parameters  X,  Y  and  Z  could  thus  be  chosen  as 


X  +  iY  =  qlq2qzql, 


4 

^  ~  ^  ^  Pnl^nl  ! 

n—1 


(lla) 

(llb) 


with  pn  being  some  constants  that  can  be  chosen  freely  to  focus  an  analysis  on  specific  properties.  Choices  of  the  pj’s 
will  be  further  investigated  shortly. 

By  taking  the  square  modulus  of  X  +  iY  and  expressing  the  l^jl^’s  in  terms  of  the  Z-coordinate  and  the  MRRs  one 
obtain  the  expression 


X^  +  Y^  =  |ging2n<73|"k4p  =  -  Zi)(Z  +  Z2)(Z  +  Z3)(Z  -  Z4), 


with  K  =  {pi  -  p2  —  P3  +  P4)  ^  and 


Zi  =  P2K3  +P3{Ki  +  K4)  -P4K1, 
Z2  =  -P1K3  +P3K2  -P4{K2  +  K4), 
Z3  =  -Pi{Ki  +  K4)  -P2K2  -P4K4, 
Z4  =  PlKi  +  P2{K2  +  K4)  +  P3K4. 


The  constants  Z^  ’s  can  be  written  in  matrix  notation  as 


/Zi 

Z2 

^3 

VZ4 


/ 

0 

K3 

KiYKi 

-Ki  \ 

(Pi 

-K3 

0 

K2 

-{K2YKi) 

P2 

-(KiYKi) 

-K2 

0 

-Ki 

Pd. 

V 

Ki 

K2YKi 

Ki 

0  / 

\P4 

(12) 


(13a) 

(13b) 

(13c) 

(13d) 


(14) 


from  which  it  is  seen  that  the  choice  of  signs  for  the  Zj's  in  Eq.  (12),  ensure  a  symmetric  matrix. 

The  coefficients  Zj  can  be  found  by  inserting  the  MRR  in  Z  leaving  only  \qj\‘^  which  can  then  be  isolated,  i.e.  for 
Zi 


Z  —  Pi\qi\'^  +P2I92P  TPsksP  +P4k4p 

=  Pi\qi\'^  +  P2{K3  -  Iqil"^)  +  P3{Ki  +  K4  -  |(7ip)  -hp4(|gip  -  Ki) 

=  {pi  -P2-P3  +P4)\qi\‘^  +P2K2  +P3{Ki  Y  Ki)  -P4K1 

=  /t-i|gi|2  +  Zi,  (15) 


while  for  Z2 

Z  =  Pi\qi\‘^  +P2I92P  +P3I93P  +P4k4p 

=  Pi{K3  -  \q2\‘^)  +  P2\q2\‘^  +  P3{-K2  +  \q2\'^)  +  P4{Ki  +  Ki  -  \q2\‘^) 

=  (pi  -P2-P3  +P4)|gip  +P1^2  -P3K2  YpiiKi  Y  Ki) 

=  -K"^|92p  -  Z2.  (16) 


By  defining 


cj,  =  X^YY‘^-  K^{Z  -  Zi)(Z  -h  Z2)(Z  -h  Z3)(Z  -  Z4)  (17) 

we  immediately  acquire  the  implicit  equation,  (p  =  0,  from  Eq.  (12),  which  yield  the  four-wave  surface  in  (X,Y,Z) 
space.  This  surface  restricts  the  evolution  of  the  fields  and  phase  to  a  two  dimensional  surface  in  a  three  dimensional 
space.  This  space  describe  the  relative  phase  of  the  fields  through  the  X  and  Y  coordinates  and  specific  field  intensities 
through  Z,  specified  by  the  choice  of  the  pj  parameters.  Choosing  pj  =  1  while  Pn^j  =  0  yield  the  evolution  of  the 
j’th  field  intensity  giving  the  exact  same  result  as  the  usual  analysis.  Choosing  instead  pi  =  P4  =  1  or  p2  =  P3  =  1 
yield  the  evolution  of  the  either  the  side  bands  or  the  inner  bands  intensities  while  choosing  pi  =  p4  =  1  =  —p2  =  —p3 
investigate  the  transfer  of  energy  between  the  inner  bands  and  the  side  bands. 
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The  Hamiltonian  in  terms  of  the  {X,  Y,Z)  coordinates  is  given  by 

ff  =  K[ri(z-  Zi)  -  r2(z  +  Z2)  -  Fsiz  +  Z3)  +  r^z  -  Z4)] 

-  y  [(Z  -  Zi)2  +  (Z  +  Z2f  +  (Z  +  Zaf  +  (Z-  Z4)2]  +  4X 

=  kAf  z  —  K(riZi  +  r2Z2  +  r3Z3  +  r4Z4) 

-  —  [4Z2  -  2(Zi  -  Z2  -  Z3  +  Z4)Z  +  Z^  +  Z^  +  Zl  +  Z^]  +  4X 
=  — Z‘^  +  K  [AF  +  i^{Zi  —  Z2  —  Z3  T  Z4)]  Z 

—  K  FiZi  +  F2Z2  +  F3Z3  +  F4Z4  +  —  (^Zf  +  Z2  +  Z3  +  Z4)  +  4A 

= -2K^Z2  +  KHZ-n  +  4A  (18) 

which  should  be  noted  to  be  independent  of  Y,  linear  in  X  and  quadratic  in  Z.  Here  the  rescaled  wave  number 
mismatch  AF  =  Fi  —  F2  —  F3  +  F4  =  ::^§r,  the  phase  mismatch  H  =  AF  +  k(Zi  —  Z2  —  Z3  +  Z4)  and  H  = 
K  [FiZi  +  F2Z2  +  F3Z3  +  F4Z4  +  ^  (Z^  +  Z|  +  Z|  +  Z|)]  has  been  defined  for  convenience.  The  dynamics  of  the 
variables  {X,  Y,  Z)  lie  on  the  orbits  formed  by  the  intersection  of  the  Hamiltonian  with  the  surface,  ^  =  0. 

The  evolution  of  the  four-wave  equations  can  of  cause  be  solved  analytically  in  the  new  coordinate  system  which 
will  be  done  in  the  following.  Any  readers  who  are  not  interested  in  rigorous  math  can  though  with  great  benefit  skip 
this  and  go  to  section  IV. 

The  Poisson  bracket  relations  for  the  new  variables  are 


(19c) 
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Using  these  the  dynamical  equations  for  the  {X,  Y,  Z)  coordinates  are 


dX 


dY 


dZ 


i{X,H} 

-2K^i{x,  z"^}  +  Kni{x,  z}  -  i{x,  n}  +  di{x,  X} 
-2K^i  {-2iZK-^Y)  +  Kfli  {-in-^Y)  -0  +  0 
-(4kZ  -  n)Y, 

-2K^i{Y,  Z‘^}  +  Kni{Y,  Z}  -  i{Y,  n}  +  Ai{Y,  X} 
-2K^i  {2iZn-^X)  +  nm  (in-^X)  -  0  +  4i  ( 

\  0 Zj 

(4«Z-f7)X  +  2||, 
oZ 

-2K^i{z,  +  Kni{z,  z}  -  i{z,  n}  +  Ai{z,  x} 

-0  +  0-0  +  4i  (in-^Y) 

-4k" 


(20a) 


(20b) 


(20c) 


In  this  form  the  coupling  between  X  and  +,  i.e.  the  evolution  of  the  phase,  consist  of  a  contribution  equivalent  to 
a  harmonic  oscillator  with  frequency  fl  and  a  power  dependent  contribution.  This  clearly  show  the  relation  that  the 
phase  of  the  four-wave  interaction  is  dependent  on  the  power. 

Differentiation  of  the  equation  for  Z,  inserting  the  equation  for  Y  and  eliminating  X  using  the  conservation  of  the 
Hamiltonian  H  gives 


d^Z 

dC2 


— 4k  ^ 


—4k  ^ 


dr 

(4kZ-H)X  +  2|| 


-k-\AkZ  -  n){H  +  2k'^Z^  -  kVLZ  +  H)  -  8k“i|| 

uZ 

-k-^[8k^z^  -  dK^nz^  +  k{ah  +  4n  +  n'^)z  -  n{H  +  n)]  -  8k“^ 

-fiz). 


d(j) 

~dZ 


(21) 

(22) 


Then,  by  using  the  common  tricks  that 


dZd^Z  d  fdZ 


dC  d^2  V 


and  that 


dZ 


giZ)  =  ^J  g{Z)dZ, 


it  is  possible  to  write  the  equation  for  Z  as  the  potential  equation 

+U{Z)  =  E, 


2  \d^ 


where 


U{z)  =  2k^Z^  -  2KnZ'^  +  2 


H  +  U- 


Z^  -  -{H  +  H)  -  8k-V(^  =  r  =  0,  Z), 

K 


(23) 


(24) 


(25) 


(26) 


and  U  is  a  constant  of  integration.  The  potential  is  thus  a  fourth  order  polynomial  and  the  potential  equation  can 
thus  be  solved  using  elliptic  equations.  The  equations  of  X  and  Y  can  thus  be  written  on  the  form 


_d  /  X 

deU 


0  -[AnZ{i)-n]\(X\  J^(Q\ 

4KZ(e)-D  0  )\Y  )^^dZ\l  ) 


(27) 
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having  a  solution  on  the  form 


X\_{  ^hom  0  \ 

Y  )  V  ^hom  ^  V  tinhorn  J  ’ 

where  Xhom7  hhom  and  Finhom  are  the  solutions  to 


d  1 

^  ^hom  ^ 

\-(  ° 

-9(0  \ 

^ ' 

-^om  J 

0  ) 

with 

diO  =  -  il, 


and 


dhinhom  _  rt  ^4* 

d^  “  ~dZ’ 


(28) 


(29) 


(30) 


(31) 


This  procedure  will  produce  the  evolution  in  the  {X,  Y,  Z)  coordinates  which  can  be  related  to  the  relative  phase 
from  the  X  and  Y  dynamics  and  the  intensities  using  the  specific  choice  of  Z  and  the  MRRs.  The  phases  of  the 
specific  field  can  then  be  found  by  splitting  the  governing  equations  for  complex  fields  Qj  into  equations  for  the  real 
amplitude  and  phases,  inserting  the  solutions  for  the  amplitudes  into  the  equations  for  the  phases  and  integrating. 

This  method  of  calculating  the  evolution  though  seem  cumbersome  and  yield  the  same  results  as  the  usual.  Another 
simple  method  giving  an  intuitive  geometric  representation  with  a  minimum  of  calculations  will  be  described  in  the 
next  section. 


IV.  THE  REDUCED  PHASE  SPACE 

In  the  last  section  the  four-wave  system  has  been  shown  to  lie  on  a  trajectory  in  the  {X,  Y,  Z)  space  and  the  method 
for  finding  the  trajectory  was  sketched.  In  this  section  another  way  of  acquiring  knowledge  of  the  dynamics  of  the 
system  is  described. 

As  shown  the  four-wave  system  was  restricted  to  the  dynamics  on  the  surface  given  hy  (j)  =  0.  This  surface  is 
determined  strictly  from  the  definitions  of  (A,  Y,  Z)  and  the  conserved  MRRs  and  not  the  physical  parameters  of  the 
system  such  as  dispersion  and  nonlinearity  and  neither  the  phases  of  the  fields.  Using  that  the  Hamiltonian  is  also 
a  constant  of  motion,  the  trajectory  of  the  four-wave  system  can  be  found  as  the  intersection  between  the  four-wave 
surface  and  the  Hamiltonian.  The  Hamiltonian  depends  on  the  dispersion  and  the  phase  of  the  fields. 

An  example  of  this  could  be  a  single  pump  parametric  amplifier,  i.e.  the  degenerate  case  at  which  q2  =  (73.  By 
choosing  pi  =  pi  =  —p2  =  —ps  =  1  the  transfer  of  energy  is  considered  since  Z  =  —1  signify  that  all  energy  is  in 
the  inner  band(s),  <72  and  <73,  while  Z  =  I  signify  that  all  energy  is  in  the  outer  bands  qi  and  (74.  As  an  example 

lets  have  the  initial  values  <71  (0)  =  <72(0)  =  qsiO)  =  and  94(0)  =  0,  i.e  for  copying  a  signal  from  <74 

to  (74  or  the  amplification  of  qi.  Such  initial  values  yield  the  four- wave  surface  shown  in  Fig.  1(a).  The  surface  is 
seen  to  be  tear-drop  shaped  between  Zmin  =  —  |  and  Z^^ax  =  1-  The  lower  bound  show  that  all  power  cannot  be 
transferred  to  the  pump  which  is  well  known  for  the  case  where  |(7i(0)|^  |(74(0)p.  The  pointy  top  of  the  tear¬ 

drop  show  the  full  conversion  of  energy  to  the  sidebands  and  the  non-analytical  shape  of  the  surface  show  that  the 
conversion  can  only  happen  in  the  infinite  time  limit.  In  the  case  where  |(7i(0)p  =  |(74(0)p  and  |(72(0)P  =  |93(0)p  the 
possibility  of  full  conversion  both  to  the  side  bands  and  to  the  inner  bands  exist,  giving  two  pointy  ends  of  the  sphere 
as  shown  in  Fig.  1(b)  with  the  initial  values  |(7i(0)p  =  |(74(0)|^  =  |  and  |(72(0)P  =  |(73(0)|^  =  |.  If  on  the  other  hand 
ki(0)P  7^  194(0)1^  and  |(72(0)|^  l<Z3(0)P  it  is  not  possible  to  have  full  power  transfer  and  thus  the  four-wave  surface 

is  smooth  as  shown  in  Fig.  1(c)  for  the  initial  values  |gi(0)p  =  |(72(0)P  =  jT,  |g3(0)p  =  ^  and  |g4(0)p  =  The 

last  example  could  of  cause  not  happen  for  the  degenerate  scheme  as,  per  definition,  <72  =  93  in  that  case. 

Lets  get  back  to  the  single  pump  parametric  amplifier.  In  order  to  determine  the  evolution  on  the  four-wave  surface 
we  only  need  to  calculate  the  Hamiltonian  and  observe  the  intersection  with  the  four-wave  surface.  For  an  amplifier 
setup  having  the  dispersion  Fi  =  0.1,  r2  =  r3  =  0.2  and  r4  =  0.3,  i.e.  having  zero  wave-number  mismatch,  AF  =  0, 
the  Hamiltonian  is  shown  in  Fig.  2(a)  and  the  intersection  between  the  Hamiltonian  and  the  four- wave  surface  is  shown 
in  Fig.  2(b).  As  seen  the  zero  wave-number  mismatch  does  not  yield  the  full  energy  transfer  as  is  well  known.  The  full 
energy  transfer  is  though  possible  for  AF  =  |(7i(0)p  -I-  |(74(0)p  which  is  seen  by  setting  Hq  =  Hopt  and  isolating  AF, 
where  the  index  0  denote  the  fields  at  ^  =  0  and  opt  denote  the  optimal  power  transfer,  i.e.  Aopt  =  0  and  Zopt  =  1- 


FIG.  2:  (a)  The  four-wave  Hamiltonian  and  (b)  the  intersection  between  the  four-wave  surface  and  the  Hamiltonian.  Both 
with  the  initial  values  gi(0)  =  ?2(0)  =  ©(O)  ~  \J \  a^nd  94(0)  =  0,  while  pi  =  P4  =  — P2  =  — Pa  =  1  and  Fi  =  0.1, 

r2  =  Fa  =  0.2  and  F4  =  0.3.  (c)  The  intersection  between  the  four- wave  surface  and  the  Hamiltonian  with  the  initial  values 
<31  (0)  =  \J 52(0)  =  53(0)  =  and  §4(0)  =  0,  while  pi  =  P4  =  — P2  =  -ps  =  1  and  Fi  =  0.1,  r2  =  Fs  =  0.2  and 
r4  =  0.5  giving  AF  =  0.2  =  |gi(0)p. 


An  example  of  full  energy  transfer  is  shown  in  Fig.  2(c)  with  the  initial  values  91  (0)  =  y  | (72(0)  =  93(0)  =  y  § 
and  54(0)  =  0,  while  Pi  =  P4  =  — P2  =  — P3  =  1  and  Fi  =  0.1,  r2  =  Fa  =  0.2  and  r4  =  0.5  giving  AF  =  0.2  =  |(7i(0)p. 
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